Static dilaton solutions and singularities in six dimensional 
warped compact ificat ion with higher derivatives 



Massimo Giovannini* 

Institute of Theoretical Physics, University of Lausanne 
BSP-1015 Dorigny, Lausanne, Switzerland 

Abstract 

Static solutions with a bulk dilaton are derived in the context of six dimen- 
sional warped compactification. In the string frame, exponentially decreasing 
warp factors are identified with critical points of the low energy /3-functions 
truncated at a given order in the string tension corrections. The stability of 
the critical points is discussed in the case of the first string tension correc- 
tion. The singularity properties of the obtained solutions are analyzed and 
illustrative numerical examples are provided. 

Preprint Number: UNIL-IPT-23 October 2000 



*Electronic address: Massimo. Giovannini@ipt.unil.ch 



1 



I. FORMULATION OF THE PROBLEM 



Consider a (4 + 2)— dimensional space-time (consistent with four dimensional Poincare 
invariance) of the form 

ds 2 = g^dxPda? = a(p)G ab dx a dx b - dp 2 - i(p)d6 2 , (1.1) 

where the (global) behavior of a(p) and 'y(p) should be determined by consistency with 
the appropriate equations of motion following, for instance, from the (4 + 2)-dimensional 
Einstein-Hilbert action This type of compactification differs from ordinary Kaluza- Klein 
schemes []],§] (see also The form of the line element given in Eq. (|1 . 1| ) can be gen- 

eralized to higher dimensions f|, although, for the present investigation, a six- dimensional 
metric will be considered. 

Most of the studies dealing with six- dimensional warped compactification assume that 
the underlying theory of gravity is of Einstein-Hilbert type. Solutions with exponentially 
decreasing warp factors have been obtained in the presence of a bulk cosmological constant 
H|7|] supplemented by either global or local |H| string-like defects (see also [[Lj] ,[n|). 



The shape of the warp factor far and away from the core of the defect is always determined, 
in the quoted examples, according to the six- dimensional Einsten-Hilbert description. Six- 
dimensional warped compactifications represent a useful example also because they are a 



non-trivial generalization of warped compactifications in five dimensions [13| . 

Up to now the compatibility of warped compactification with theories different from 
the Einstein-Hilbert description have been analyzed mainly in the five-dimensional case. In 
particular the attention has been payed to gravity theories with higher derivatives |I4]-[T6|. 
In [pjj the interesting problem of the compatibility of a five-dimensional compactification 



x The conventions of the present paper are the following : the signature of the metric is mostly 
minus, Latin indices run over the (3 + l)-dimensional space, whereas Greek indices run over the 
whole (4 + 2)-dimensional space. 
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scheme with higher derivative gravity and dilaton field has been investigated. In |18j the 
simultaneous presence of higher dimensional hedgehogs and of higher derivatives gravity 
theories has been discussed in a seven-dimensional space-time. 

The simultaneous presence of dilaton field and higher order curvature corrections is rel- 
evant for different reasons. There are two possible approaches to warped compactification 
schemes. In the first approach it is assumed, from the very beginning, that warped compact- 
ifications have nothing to do with string theory. The problem then reduces to the search 
of a suitable field theoretical model of higher dimensional (global or local) defects living 
along the transverse dimensions. A complementary approach is to postulate that warped 
compactifications are (somehow) connected to string theory. Even in this second approach, 
however, the dilaton field and the (possible) higher curvature terms (arising from the string 



tension corrections) are usually neglected [TI|. The assumption that the dilaton field is 



strictly constant and that the string tension corrections are vanishing should be relaxed in 
order to investigate the (possible) interplay of warped compactification schemes with string 
models. 

Five-dimensional warped compactifications have been studied in the context of general- 
ized gravity theories with higher derivatives [P^IB|. The inclusion of the dilaton field poses, 
however some problem. The main results concerning the interplay between five-dimensional 
warped compactifications and string motivated effective action can be summarized as fol- 
lows jn|]. If the dilaton is constant (but higher order curvature corrections are still present) 



five-dimensional warped compactifications are consistent with the effective action derived 
from the string amplitudes corrected to first order in the string tension. 

If the dilaton is linear in the bulk coordinate warped compactification is not consistent 
with the presence of string tension corrections unless the conformal coupling between the 
dilaton and the string tension corrections is absent [] |L7| . Interesting solutions interpolating 



2 The effective action used in order to derive these results has been studied Einstein frame. 
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between a naked singularity and a warped regime (for large bulk coordinate) have also been 
obtained. The purpose of the present investigation is to is to analyze the solutions of the 
metric ( |1.1| ) in the context of the string effective action [0,0 with and without higher 



derivatives corrections fElLE2l (see also [23,2 



Up to now the compatibility of six- dimensional warped compact ification with gravity 
theories inspired by string amplitudes has not been analyzed. In six- dimensional warped 
compactifications the transverse space is not flat (in contrast with the five-dimensional 
case). This observation implies on one hand that the internal space is larger and, on the 
other hand, that two warp factors are generically allowed. The tree-level solutions have more 
general singularity properties if compared to the five- dimensional case. These solutions are 
"Kasner-like" and they are the static analog of their time-dependent counterpart which is 



often discussed in cosmological solutions [25]. The inclusion of string tension corrections 



introduces also differences with respect to the five- dimensional case. 

In six dimensions exponentially decreasing warp factors correspond to critical points of 
the (3 functions computed at a given order in the string tension. Define, in fact, 7i = d p (In a) 
and T = <9 p (hi7). By critical points we mean those (stable or unstable) solutions for 
which 7i and T are simultaneously constant and negative. The critical points of the system 
correspond to a static dilaton fields which either increases or decreases (linearly) for large 
p. Also the constant dilaton solution is possible. This situation should be contrasted with 
the five-dimensional case where only one warp factor is present. 

The critical points are not always stable. If the initial conditions of the dilaton and of 
the warp factors are given around a given critical point (say for p = p ) it can happen that 
for larger p the compatibility with the /3-functions will drive the solution away from the 
(original) critical point. In this case singularity may also be developed. 

The present analysis is not meant to be exhaustive and suffers of two obvious limitations. 
In order to make an explicit calculation the effects of the first string tension correction has 
been studied. This is just an example since, when singularities are developed, all the string 
tension corrections should be included. The second point to be emphasized is that possible 



corrections in the dilaton coupling have also been neglected. This might be justified in 
some regimes of the solutions but it is not justified in more general terms. With these two 
warnings in mind the reported results should be understood more as possible indications 
than as a firm conclusion. 

The present analysis has been performed in the string frame where the dilaton field is 
directly coupled both to the Einstein-Hilbert term and to the first string tension corrections. 
It is interesting to notice that singularity properties of a given solution may change from one 
frame to the other. A linearly increasing dilaton (in the string frame) results in curvature 
singularities in the Einstein frame (and vice- versa). The Einstein and the string frame are, 
however, physically equivalent. Suppose that the curvature invariants are regular in one 
frame and suppose that the dilaton is a smooth function of the bulk radius (interpolating, 
for instance, between two constant regimes). Then, the geometry looks regular in both 
frames enforcing the physical equivalence of the two descriptions. 

In order to assess that six-dimensional warped compactification is fully consistent with 
the presence of the dilaton field and of string tension corrections some requirements should 
be, in our opinion, satisfied. The goal would be to obtain stable warped solutions with 
constant dilaton. Stable warped solutions means that the critical point of the low energy f3 
functions (truncated to first order in the a' corrections) lead to exponentially decreasing warp 
factors. The constancy of the dilaton field implies the constancy of the coupling constant 
(i.e. exp [0/2]) and this is a necessary (even though not sufficient) condition for the stability 
of the solutions. Notice that once the dilaton field and the string tension corrections are 
simultaneously present the problem of the singularity of the geometry cannot be disentangled 
from the problem of the dilaton relaxation. 

The plan of the present investigation is then the following. In Section II the tree-level 
solution of the low energy string effective action will be studied. In Section III the first 
string tension correction will be included and the corresponding equations of motion will be 
derived. In Section IV the critical points of the obtained dynamical system will be analyzed. 
In Section V the stability of the critical points will be scrutinized and some numerical 
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examples will be presented. Section VI contains some concluding remarks. In the Appendix 
useful technical results are reported. 



II. TREE-LEVEL SOLUTIONS 

If the curvature of the background is sufficiently small in units of the string length 3 
(denoted by A s ) the massless modes of the string are weakly coupled and the dynamics can 
be described by the string effective action in six dimensions |T9|j2( 



= j d 6 xV=g£ {0) = -2^4 J d 6 xV=ge-* [r + g^datdpf + A] , (2.1) 

where is the dilaton field, X s = \fa' (a' being the string tension). Eq. ( |2.1| ) is written 
in the string frame where the string scale is constant and the Planck scale depends upon 
the value of the dilaton coupling, i.e. g(<p) = e^ 2 ! . If the dilaton coupling is constant the 
the string frame coincides with Einstein frame. If, as in the present case, = 4>(p) the two 
frames are equivalent up to a conformal transformation. In the present paper the string 



frame will be used. In Eq. (2Jj) the minimal field content (i.e. graviton and dilaton) has 
been assumed together with a bulk cosmological constant A. The effective action (|2.1| ) is 
derived by requiring that the usual string scattering amplitudes are correctly reproduced 
to the lowest order in a'. The requirement that the equations of motion derived from Eq. 
( P-1P are satisfied in the metric ( |1 . 1| ) is equivalent to the requirement that the background is 



conformally invariant to the lowest order in a'. Notice that in eq. (|2.1 ) the contribution of 
the antisymmetric tensor field has been neglected. This is justified within the spirit of the 
present analysis but this might not be justified in more general terms. We will come back 
on this point in the following Sections. 



The equations of motion derived from the action of Eq. (|2.1| ) can be written as 

R - g a(3 d a <pd^ + 2g^V a Vp(j> + A = 0, (2.2) 



3 In the discussion string units (i.e. X s = 1) will be often used. 
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Rp* - ^g^R + ^g^g aP 9 c 



A 

= °- 



(2.3) 



By now using the metric of Eq. (|1.1Q the explicit form of the (a, a) and (0, 9) components 
of Eq. Q will be Q 



2A + J ^^2 + 3J ^^ 7^ + 6 7^ 2 + 2J ^^ ' + 6 7^ , -2.F0'-6ft0' + 2 /2 -4 0" = 
2A + Km 2 + 8ft'-8ft0' + 20' 2 - 40" = 0, 



(2.4) 
(2.5) 



whereas the explicit form of the of Eq. (|2.2j) will be 



2A- T 2 WH 2 -2T' - %H' + 2F<1>' + 8H(p' -20' 2 + 40" = 0. (2.6) 



Notice that Eq. (|2.2j ) has been multiplied by a factor two and Eq. fl2.3|) has been multiplied 
by a factor of four in order to get rid of rational coefficients. In deriving Eqs. (|2.4|)-([2.6|) 



it has been assumed that, in Eq. (|1.1| ) G a b = rj ab where r] ab is the Minkowski metric. Eqs. 



^-(lA)) admit exact solutions whose explicit form can be written as: 

2a 



<r(p) = tanh 

7(p) = tanh 
e m = _| siim 



-A 



2 



(p - Po) 
(P - Po) 



20 



v^A 



(P - Po) 



4q+/5-1 



cosh 



V ^A -.,-4a-/?-l 

— S~ (P-Po) 



(2.7) 
(2.8) 
(2.9) 



The exponents a and (3 satisfy the condition 4a 2 + (3 2 = 1. Eqs. (|2~7j ) -(|2T9| ) lead to a 



physical singularity for p — > po- All curvature invariants associated with the metric ( |1 . 1| ) for 



the specific solution given in Eqs. ( |2.7] )- fl2~9| ) diverge for p — > po (see Appendix A for the 
details). If a = f3 the Weyl invariant vanishes but the other invariants are still singular. 
If A -> 0, Eqs. M-(pl) are solved by 



4 As previously mentioned, for convenience, the following notations are used TL = (In a)', T 
(In 7)' where ' = d p . 
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a(p) = (p - p ) 2a , (2.10) 

7(P) = (P ~ Po) 2/3 , (2.11) 
= (p - p )- 1+ia+ ^ (2.12) 

provided Aa 2 + (3 2 = 1. Eqs. (p77|)-(F^D and (plUH P^) are Kasner-like Q solutions 
whose time-dependent analog has been widely exploited in the context of string cosmological 
solutions |25| . 



In the case of A < this system of equation has a further non trivial solution which 
is given by 7i, T and <$' all constant. In this case the solution of the previous system of 
equations is given by 

4A + 4H 2 + ^ 2 = 0, 2(f)' = AH + J 7 . (2.13) 



In the particular case where a(p) oc 7(p) ( i.e. H = J 7 ), the solution is H = — y— 4A/5. As 
we will discuss in Section IV this solution is not always stable. 

The presence of curvature singularities in the tree-level solutions obtained in Eqs. ( |2.7| )- 
( p.9| ) and ( |2.10| )-( f2.11| ) suggests that there are physical regimes where the curvature of the 



geometry will approach the string curvature scale. In this situation higher order (curvature) 
corrections may play a role in stabilizing the solution and should be considered. The follow- 
ing part of this investigation will then deal with the inclusion of the first a' correction. This 
analysis is of course not conclusive per se since also higher orders in a' should be considered 
as it has been argued in Section I. 



5 For truly Kasner solutions the sum of the exponents (and of their squares) has to equal one. 
In the present case only the sum of the squares is constrained and this is the reason why these 
solutions are often named "Kasner-like". 
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III. FIRST ORDER a' CORRECTIONS 



Consider now the first a' correction to the action presented in Eq. Q2.1|) . The full 
action is, in this case [0] 

1 



S = S<® + 



2\t 



d x v '—g e 



R + 9 af) d a <i>dp<t> + A - e R^ aP Rr«P , (3.1) 



where e = ka'/A and the constant k takes different values depending upon the specific theory 
( k = 1 for the bosonic theory, k = 1/2 for the heterotic theory). Notice that the assumption 
made in the previous Section (concerning the absence of antisymmetric tensor field) reflects 
in the fact that the first a' correction appears precisely in the form reported in eq. ( |3.1|) . If 
the antisymmetric tensor field would not be vanishing the first string tension correction to 



the tree-level action would look like [EI 



rA 

fjtv 



, ( E> Tfiwaf) 1 Tyavaf) tt tt> 



(3.2) 



where H^a is the antisymmetric tensor field strength. If vanishes identically the 

only correction to the tree-level action is the one coming from the first term of Eq. Q3.2|) . 
However, if H^ ua ^ 0, the situation changes qualitatively. We leave this intriguing issue for 
future investigations. 

The fields appearing in the action (|3.1|) can be redefined (preserving the perturbative 
string amplitudes) pTjpBI . In order to discuss actual solutions it is useful to perform a 
field redefinition (keeping the a model parameterization of the action) that eliminates terms 



with higher than second derivatives from the effective equations of motion [21,22] (see also 
p4]j). In six dimensions the field redefinition can be written as 



+ e 



R+(3 + 2n)g^g^d a <f ) d p <\> 



(3.3) 



where n (the number of transverse dimensions) is equal to 2 in the case of the present 
analysis. Dropping the bar in the redefined fields the action reads : 
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s 



2A! 



(Fx, 



-ge 



R + g^d^d^ + A - e[Rl GB - g^g^d^drfdrfd^] , (3.4) 



where 



p2 p r>iivaf3 a p T3V- V i E>2 

-^EGB _ tiiivaflti — QK^K + -ft , 



(3.5) 



is the Euler- Gauss- Bonnet invariant (which coincides, in four space-time dimensions, with 
the Euler invariant Q). The equations of motion can be easily derived by varying the action 
( p.4| ) with respect to the metric and with respect to the dilaton field (see Appendix B for 
details). The result is that 

2 A + F 2 + 3 FH + GH 2 + 2 F' + QH' -2T$ - 6ft0' + 2 0' 2 - 40" 
+e [-3 T 2 H 2 -9FH 3 - 3ft 4 - 6ft 2 T' - 12FHH' -QH 2 H'+ 
§T 2 H<\> + 2AFH 2 0' + 18ft 3 0' + Y2HT' 0' 
+ 12^ft'0' + 24ftft'0'- 12^ft0' 2 - 12ft 2 0' 2 + 



2 itn 



2 0' 4 + 12.7^0" + 12 ft 
2A + 10ft 2 + 8ft'-8ft0' + 2 0' 2 -40" + 
e [- 15 ft 4 - 24 ft 2 ft' + 48 ft 3 0' + 48 ft ft' 0' - 24 ft 2 0' 2 + 



+ 24 ft 



2 



-2A-J zr2 -4J !r ft-lOft 2 -2.r-8ft' + 2j !r 0' + 8ft0'-2 0' 2 + 4. 
+e [6^ 2 ft 2 + 24J^ft 3 + 15ft 4 + 12ft 2 .F' + 24.7 r ftft' + 24ft 2 ft' 
-4 T 0' 3 - 16 ft 0' 3 + 6 0' 4 - 24 0' 2 0"1 = 0. 



(3.6) 



(3.7) 



(3.8) 



In the limit e — > the equations derived in Section II are recovered. Eqs. ( |3 . 6| ) — ( pT8|) can 
be studied in order to investigate two separate issues. The first one is the existence of 
critical points leading to exponentially decreasing warp factors. The second issue would be 



6 The Euler-Gauss-Bonnet invariant is particularly useful in order to parameterize quadratic cor- 
rections in higher dimensional cosmological models [p6-28]. 
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to analyze the stability of these critical points. Indeed, in the solutions given in Eqs. 
( p.9| ) a singularity is developed. This result is based, however, on the tree-level action. If 
the first a' correction is included this feature might change. Moreover, if A = there are no 
warped solutions at tree-level. These solutions might emerge, however, when a' corrections 
are included. For instance, the "Kasner-like" branch of the solution might be analytically 
connected to a "warped" regime where Ti and T are both constant and negative. These 
questions will be the subject of the following Section. 

In closing this section we want to recall that in [29f it has been claimed that (five- 



dimensional) warped compactifications cannot be obtained in the low energy supergravities 

coming from string theory. In this context it has also been suggested that higher order 

corrections to the Einstenian action might help in evading this conclusion. In the following 
we will explore an analogous possibility. 

IV. DYNAMICAL SYSTEM AND CRITICAL POINTS 

Defining 

x(p)=H{p), y{p)=T{p), z{ P ) = <t>'{p), (4.1) 
Eqs. (|3.6|) - (|3.8D can be written as 

Px(j>) X' + Py{p) V' + Pz{p) Z* + Po(p) = 0, (4.2) 

q x { P )x , + q z (p)z' + q o {p) = 0, (4.3) 

w x (p) x' + w y (p) y' + w z (p) z' + w (p) = 0. (4.4) 

The p are 

p x {p) = 6{1 - e[x(x + 2y) - 2z{y + 2x)}}, 

p y (p) = 2[1 - 3ex{x - 2z)}, p z {p) = -4[1 - 3ex{x + y)}, 

p (p) = 2A + y(y + 3x) + 6x 2 - 2z{y + 3x) + 2z 2 

+e[-3x 2 (y 2 + 3xy + x 2 ) + 6x(x + y)(y + 3x)z - I2x{x + y)z 2 + 2z\ (4.5) 
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whereas the q are 

q x {p) = 8[1 - 3ex(x - 2z)}, q z {p) = -4(1 - 6ex 2 ), 

q (p) = 2A + 2(z 2 - Axz + 5x 2 ) + e[-15x 4 + A8x 3 z - 2Ax 2 z 2 + 2z 4 ], (4.6) 

and finally the w are 

^x(p) = -8[1 - 3tx(x + y)], 

w y (p) = -2[1 - Qex 2 }, w z (p) = 4(1 - 6ez 2 ) } 

w (p) = -2A - (y 2 + Axy + 10x 2 ) + 2(y + Ax)z - 2z 2 

+e[3x 2 (2y 2 + 8xy + 5x 2 ) - 4(y + 4x)^ 3 + 6z 4 ]. (4.7) 

Eqs. (4.2)— (]4^) can be analyzed both analytically and numerically. Interesting analytical 



conclusions can be obtained, for instance, in the case x(p) = y(p)- In the remaining part of 
the present Section the stability of the system will be analyzed. Physical (i.e. curvature) 
singularities will be investigated. The logic will be, in short, the following. The critical 
points will be firstly studied analytically. Indeed, it happens that, in spite of the apparent 
complications of the system, there are limits in which analytical expressions can be given. 
Numerical examples will be given for the other cases which are not solvable analytically. 
The examples studied in the present Session will be exploited in Section V for some explicit 
numerical solutions. 



The critical points of the system derived in Eqs. (|4.2j )- (|4.4j ) are defined to be the one 
for which [3(J 

x'(p)=0, y'(p)=0, z'(p) = 0. (4.8) 

The critical points interesting for warped compactifications are the ones for which x(p) and 
y(p) are not only constant but also negative. These points lead, for large p, to exponentially 
decreasing warp factors in the six- dimensional metric. 

In the case of Eq. ( |4.8| ) the solution of the system reduces to a system of three algebraic 
equations in the unknowns (x, y, z). The three equations of the system are simply given by 
the conditions 
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Po(p) = °) Qo(p) = °> w o{p) = 0, 



(4.9) 



namely, 



2A + y(y + 3x) + Qx 2 - 2z(y + 3x) + 2^ 2 

+e[— 3x 2 (y 2 + 3xy + x 2 ) + 6x(x + y)(y + 3x)z — 12x(x + y)z 2 + 2z 4 } 
2A + 2(z 2 - 4xz + 5x 2 ) + e[-15x 4 + 48x 3 z - 24x 2 z 2 + 2z A ] = 0, 
-2A - {y 2 + Axy + 10x 2 ) + 2{y + Ax)z - 2^ 2 
+e[3x 2 (2y 2 + 8xy + 5x 2 ) - 4(y + 4x)z 3 + 6z*} = 0. 



(4.10) 
(4.11) 

(4.12) 



By summing up, respectively, Eqs. ( 4.10j ) and ( |4.11| ) with Eq. ( 4.12|j the following two 
algebraic conditions are obtained: 

(4x + y - 2z)[3ex 3 - Aez 3 + 3ex 2 (y + 2z) + x{-l + deyz)} = 0, 

(Ax + y - 2z)[y(-l + 6ex 2 ) - Aez(z 2 - 3x 2 )] = 0. (4.13) 

A third algebraic relations is obtained by subtracting Eq. (|4.11|) from Eq. (|4.10|) : 



(x - y)(Ax + y- 2z)[-l + 3ex(x - 2z)\ = 0. 



(4.14) 



If x = y Eqs. ( f4.10|) -( f4.12|) are solved provided 

Ax + y 

x = y, z = — - — , 
16A + 20x 2 + 265ex 4 = 0, 



(4.15) 
(4.16) 



and an explicit solution corresponding to the critical points of the system ( |4.10|) -( [i.l2 ) can 
then be written as 



<7(p)=eW, 7 (p)= 70 (T(p) 
z pi 

where 7o and O are integration constants and where 



(4.17) 



(4.18) 
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is obtained from the real roots of the algebraic equation appearing in Eq. ( f4.16|) . A priori, 
depending upon the sign of A and upon the relative weight of Ae, pi can be either positive 
or negative. 

The solution derived in Eqs. ( |4 . 1 6| ) does not exhaust the possible critical points even in 
the case x = y. Consider, in fact, Eqs. ( |4.1U| )-( |4.12| ) in the case x = y but with z ^ 5x/2. 
In this case Eq. (|4.10| ) and Eq. ( [4.1 Ij ) lead to the same algebraic equation 



2A + 2 (5x 2 -Axz + z 2 ) +e (-15x 4 + 48x 3 z - 24x 2 z 2 + 2 z A ) = 0, 



(4.19) 



whereas Eq. (|4.12j ) leads to 



2A - 15x 2 + 10xz -2z 2 + e (45x 4 - 20xz 3 + Q z A ) = 0. 



(4.20) 



Eqs. ( |4.19| )- (|4.20| ) cannot be further simplified and the best one can do is to solve them 
once the values of e and A are specified. Ifj] A = —1 and e = 0.1 the real roots of the system 
of Eqs. (^T9|)-(^20|) are given by: 



(xi, z x ) = ±(0.6974, 1.7436), {x 2 , z 2 ) = (±1.3775, ^2.3429), 
(x 3 , z 3 ) = (±0.1338, =f0.7176), (x 4 , z A ) = (±1.9645, ^0.5744). 



(4.21) 
(4.22) 



On top of the real roots reported in Eqs. ([4.21|) - ([4.22|) there are also imaginary roots which 



are not relevant for the present discussion. The first pair of roots in Eq. ([4.2 1|) can also 
be obtained from Eq. ( |4.16|) . Notice, also, that {x\, z±) is the only root where the z and x 
can be simultaneously negative (when the minus sign is chosen). The other roots are such 
that x and z always have opposite signs. The same discussion can be repeated for different 
values of A and e. For instance, if A = — 2 and e = 0.1 the real roots of Eqs. ( ^4.19|) - ([4.20| ) 



arc 



7 The numerical roots reported in the present session are not just illustrative. They will be 
exploited in the numerical solutions discussed in the following Sessions. 
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(xi, zi) = ±(0.8857, 2.2143), (x 2 ,z 2 ) = (±1.3490, ^2.3077), (4.23) 
(x 3 , z 3 ) = ±(1.0275, 1.6067), (ar 4 , z A ) = (±0.2455, ^0.9111), (4.24) 
(x 5 , x 5 ) = ±(0.8563,0.7023), (x 6 , z Q ) = (±1.9285, ^0.5468). (4.25) 



If x ^ y and, simultaneously, z ^ (Ax + y)/2 the system of Eqs. d4.1U| )— ( fi. 12| ) can still 



be solved. From Eqs. ( [4. 13|) and Q4.14|) an ansatz for the algebraic solution can be written 

as 

3cx 2 - 1 



z 



6 ex 

1 - 9ex 2 - 81eV ± 297e 3 x 6 
54e 2 x 3 - 324e 3 x J 



By inserting this ansatz back into Eqs. ( |4.10| )— fl4.12|) the same equation for x is obtained 
from all three equations of the system, namely, 

5 1 1 7x 2 25 ex 4 n , A nr ^ 

2A+ nl + mJT* + W^ + — + — = - (427) 



Thus, the critical points of the system of Eqs. ( |4.10| )-( |4.12| ) [in the case x ^ y, x ^ and 
z 7^ (Ax + y)/2] are given by 

_ 3ex 2 - 1 _ 1 - 9ex 2 - 81eV + 297e 3 x 6 
Z ~ 6ex ' y ~ 54e 2 x 3 - 324e 3 x 5 ' ( ' ^ 

. . 1 x 2 ( , 5 \ , 7x 6 25 ex 8 , , 

P x,A,e = -± -± 2A+ x 4 ± + = 0. 4.29 

v ' ' ; 648 e 3 27e 2 \ I2t) 3 8 V ; 

Eq. (|4.29| ) should be solved for real values of x. The critical points corresponding to negative 
(and real) roots of Eq. ( fl.29| ) lead to exponentially decaying cr(p). Once the real roots of 
Eqs. ( |4.29| ) are obtained, Eqs. ( [4.28| ) will give the wanted values of y and z. Eq. ( [4.29| ) 



does not always have real roots. In Fig. [IJ P(x, A, e) is reported as a function of x and 
for different values of A and e. The specific values are only representative and will be used 
in specific numerical solutions later on. For instance, it can be argued from Fig. [l| that if 
e = 0.1 real roots of Eq. (|4.29|) start appearing for A < —6. 
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FIG. 1. P(x,A, e) appearing in Eq. Q4.29|) is plotted in the case of e = 0.1 (left) and e = 1 
(right) and for different values of A. In the left plot the thick line corresponds to the case A = —6 
where P(x, A, e) has real roots. The other thin lines (from bottom to top) correspond, respectively, 
to A = —5, —4, —3, —1. In the right plot the thick curve correspond to the case A = —2 and the 
two thin lines correspond, from bottom to to top, to A = —3 and A = —1. 

V. STABILITY OF THE CRITICAL POINTS, SINGULARITIES AND 

NUMERICAL EXAMPLES 



If x(p) = y(p) [i.e. 7(p) = jocr(p)], Eqs. (|4.2| )-(fOD lead exactly to the same equation so 
that the only two independent equations obtained in this case are: 



where 



and 



a x (p) x' + a z (p) z + a (p) = 0, 
b x (p) x' + b z {p) z' + b (p) = 0. 



a x (p) = 8[1 - 3ex(x - 2z)], a z (p) = -4(1 - 6ex 2 ), 
b x (p) = -10(1 - 6ex 2 ), b z {p) = 4(1 - 6ez 2 ), 



ao(p) = [2A + 10x 2 - 8xz + 2z 2 + e(-15x 4 + A8x 3 z - 2Ax 2 z 2 + 2z% 
b (p) = [-2A - 15x 2 + lOxz - 2z 2 + e(45x 4 - 20xz 3 + 6z% 



(5.1) 



(5.2) 



(5.3) 
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Notice that the If A = 0, from Eq. ( 4.16|) p\ becomes purely imaginary. If A < and if p\ 



has to be real there are two relevant critical points, namely 



"'=^/8i:V 1+ V 1+ f A - £ ' (5 - 4) 

where, for convenience, A = — A_ (with A_ > 0). In order to have a warped compactification 
scheme, pi < and, therefore, only one critical point satisfies this requirement. 

If, as a warm-up, e^Oa critical point of the system is indeed given by Eq. (|2.13| ). The 



tree-level critical point is not stable. Indeed, in the case e = Eqs. ( |5.1|) can be written as 

x' = f x (x, z), z = f z (x, z), (5.5) 

where 

5 

fx\X, z) xz —x , 

h(x,z) = ^-^x 2 + j. (5.6) 

The partial derivatives of f x (x,z) and f z (x,z) (i.e. d x f x , d z f x and d x f z , d z f z ) define a 
matrix whose entries should be evaluated in the critical points, corresponding, in the present 
case, to z = (5/2)x, x = — a/4A__/5. The eigenvalues of this matrix are simply given by 
<5i ; 2 = =Fa/5A_/4. Therefore, the tree-level critical point is an unstable node since the 
eigenvalues are both real and with opposite signs (i.e. <5i<5 2 < 0). 

Suppose now to turn on the first a' correction. Then e ^ 0. In this case the tree-level 
discussion can be repeated. If e 7^ the dynamical system will now become [always in the 
case x(p) = y(p)} 

x' = g x (x,z), z'=g z (x,z), (5.7) 



where now, from eq. ( p7T 



a z b - a b 

9x{x,z) 



b x a z — a x b 



z 

1 



X 



9z (x,z) = f°-_ aob * (5.8) 
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In analogy with the previous case the matrix of the derivatives can be constructed. The two 
eigenvalues of such a matrix are complex conjugate numbers of the form 



W h2 = VA-Hk)±i v(k)]. (5.9) 

where k = A_e and the subscripts 1 and 2 correspond, respectively, to the plus and minus 
signs. The expressions are quite cumbersome so that they are reported in the Appendix C. 
In order to determine the stability of the system the sign of n(k) is crucial. It turns out, 
from a numerical analysis of Eqs. (|C.1|) -( |C^ ), that 

H(k) > 0, for k = A_e > 0.65, 

fi(k) < 0, for k = A_e < 0.65. (5.10) 

Thus, for A_e > 0.65 there is an unstable spiral point whereas for A_e < 0.65 there is a 
stable spiral point [[3(] . 



In closing this session an interesting solution can be mentioned. In the case A = and 
z = (constant dilaton case) the relevant equations can be written as : 

6[1 - ex(x + 2y)]x' + 2[1 - 3ex 2 }y' + 2A + 6x 2 

+y{y + 3x) - 3ex 2 {y 2 + 3xy + x 2 ) = 0, (5.11) 
8(1 - 3ex 2 )x' + 2A + 10x 2 - 15ex 4 = 0. (5.12) 

The critical point of the system can be obtained also in this case and it corresponds to x c = 



y c = — a/2/(3 e). The system can then be written in its canonical form, namely, x' = F x (x, y) 
and y' = F y (x, y). The eigenvalues of the matrix of the derivatives evaluated in (x c , y c ) are, 
respectively, 5i = 5/\/6e and 5 2 = 5(3e — l)/\^Qe. Notice that 5i5 2 = [25(3e — l)]/(6e). This 
means that the critical point is a saddle point for e < 1/3 and an unstable node for e > 1/3 
p0| . This solution, even though interesting, is only illustrative : the presence of the inverse 
coupling in the critical point clearly points towards non-perturbative effects which should 
be properly discussed by adding higher orders in a'. 
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8 A =0.1 



- 1.7436' 




FIG. 2. A numerical solution of the system given in Eq. ( |5.1| ) is reported in the case A = — 1 
and e = 0.1. The solution interpolates between the Kasner-like regime (for p — ► 0) and the critical 
point given by the solution of Eq. (|4.16| ) which is explicitly given (for the chosen values of A and 
e) by the first pair of roots appearing in Eq. ( fl.21 ). 



A. Numerical examples 



A Kasner-like branch can be analytically connected to a constant curvature solution 
where x < and the dilaton is linearly decreasing (z < 0). Suppose that A = — 1 and 
e = 0.1. Then, we can see that there exist a solution connecting the Kasner-like regime to 
a critical point given by Eqs. (|4.16|) and (|4.21| ). The numerical result is reported in Fig. 
0. For large p the solution matches exactly the critical point which can be obtained (with 
A = — 1 and e = 0.1) from Eq. (|4.16|) and which are reported in the first pair of roots in Eq. 

If the initial conditions are not given in the Kasner-like regime but in the vicinity of 
a critical point there are various possibilities. It can happen that the system is driven 
(smoothly) towards another critical point or it can happen that a singularity is reached. 

Suppose that the initial conditions for the numerical integration of the system of Eq. 
( pM| ) are given in in p ~ requiring that x(0) < and z(0) < according to the solution 
of Eqs. ( [4.161 ). In this case, the system seats exactly on a critical point whose properties 
have been previously analyzed. In the case A = — 1 and e = 0.1 the initial conditions are 
fixed for x(0) = X\ and z(0) = z\ where {x\ y z\) are given by Eq. Q4.21| ). In this case 
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FIG. 3. The numerical integration of Eqs. ( p.l| ) is reported in the case A = —1 and e = 0.1. 



The critical points in this specific case were given explicitly in Eqs. ( 4.21 )— ( 4.22 ). Initial conditions 
are given in one of them, namely (xi, z\). 

the system evolves towards a different critical point, namely it happens that x(oo) — > x% 
and z(oo) — > z$ where (x^, z$) are given by Eq. (|4.22|) . This conclusion can be obtained 
by direct numerical integration. In Fig. ^ the numerical integration is reported for the 
mentioned initial conditions. It can be seen that x(50) = 0.133 and z(50) = —0.7176 which, 
indeed, coincide with x^ and z% within the numerical accuracy of the present analysis. The 
singularity structure of the curvature invariants can be appreciated by analyzing the behavior 
of 



h(p) = R^RT *, Hp) = C, uaP C^ 
h(p)=R, h{p) = R^RT 



(5.13) 



whose general form is reported, for the metric (1.1), in Appendix C. Notice that if x = y, 
as in the example of Fig. ^, the Weyl invariant vanishes. 

Suppose now that A = —2 and e = 0.1. The critical points are given by Eqs. (|4.23 )— 
( |4.25| ). The initial conditions imposed in p = will be such that x(0) < and z(0) < as 
in the previous numerical example. Requiring, from Eq. (|4.23|) , x(0) = X\ and z(0) = Z\ the 
numerical integration leads to a divergence for large p. In this specific case, the singularity 
occurs around p ~ 6. Again, the proof of this statement can be given by direct numerical 
integration whose results are reported in Fig. [|. Numerical examples concerning the case 
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FIG. 4. In the case of A = —1 and e = 0.1 the curvature invariants corresponding to the 
solution reported in Fig. || are reported. The Riemann invariant appears in the left plot with full 
thick line. At the right the Ricci invariant is reported with dot-dashed line. The scalar curvature 
is reported (always in the right plot) with the full (thick) line. 
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FIG. 5. In the case A = —2 and e = 0.1 the numerical integration of the system ( |5.1| ) is 
reported. The initial conditions are given in a critical point of the system , namely the point 
(xi, z\ ) reported in Eq. ( |1.23 ). 
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FIG. 6. The Riemann (left) and Ricci (right) invriants are reported for the numerical solution 
obtained in Fig. |5| 

x 7^ y will now be given. The specific values of e and A are just meant to illustrate some 
interesting features of the solutions. As a general remark we could say that in the case x 7^ y 
singularities seem to be more likely. This statement is only based on a scan of the numerical 
solutions for different values of A and e. 

Suppose that e = 0.1. Thus, from the examples of Fig. (P, A needs to be sufficiently 
negative in order to produce real (negative) roots in Eq. ([4.29|) . Suppose, for instance that 
A = -10. Then, from Eqs. (|]2|)-(j42|), the critical points will be, 



(x u y x , z x ) = (-2.0258, 0.6386, -0.1902), 
(a?2, 2/2, z 2 ) = (-0.6889, -3.3432, 2.0748). 



(5.14) 



If the initial conditions for the warp factors and for the dilaton are set in (x\, yx, zx) of 
Eq. ( |4.25| ) then the results of the numerical integration are reported in Fig. [7]. A physical 
singularity is developed at a finite value of p where both the Weyl and the Riemann invariant 
blow up. The same kind of solutions can be obtained if the initial conditions are set not in 
(x 2 , y2, z 2 ) of Eq. ( [5.141 ). Also in this second case a singularity is developed at finite p. 

In order to describe the qualitative behavior of the system the value of e can be changed. 
From Fig. |l] (left plot) A = —2 can be selected and from Eqs. ( 4.28 )-( 4.29 ) the critical 
points are found to be 

(xx, yx, zx) = (-0.8679, 0.6052, -0.2419), 
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FIG. 7. In The case of A = —10 and e = 0.1 the numerical integration of the system ( |4.4| ) is 
reported in the left plot. The evolution of x (dashed line), y (full thick line) and z (thin line) are 
illustrated when the initial conditions are given in (x\,y\,Z\) of Eq. ( [4.25|) . In the plot at the right 
the Weyl (dashed line) and the Riemann (full line) are reported. 

O2, 2/2, z 2 ) = (-0.1638, -3.5390, 0.9350). (5.15) 

In this case the results of the numerical integration are reported in Fig. [S[ Also in this 
example a physical singularity is developed at finite value of p. 

In the present Section we examined systematically the critical points of the dynamical 
system arising from the low energy f3 functions truncated to first order in the string tension 
corrections. For practical purposes two cases have been distinguished namely the one where 
the warp factors are proportional and the one where the warp factors are both exponentially 
decreasing but not proportional to each other. From the analysis of both cases we can 
conclude that there are no (stable) critical points of the /3-functions where the dilaton is 
frozen. This means that stable (or unstable) critical points with decreasing dilaton field 
cannot be analytically connected with solutions where the dilaton coupling is frozen to a 
constant value. However, a large class of critical points leads to warped compactifications 
in the sense that the warp factors are exponentially decreasing and the dilaton coupling is 
also exponentially decreasing (without getting fixed to a constant value). 

As pointed out in the introduction, the present investigation suffers from various limita- 
tions. That is why the present results can only be considered as illustrative. However, it is 
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FIG. 8. In the left plot the numerical integration of eqs. ( |4.4[ ) is illustrated for A = —2 and 
e = 0.1. The notations are exactly the same as in Fig. [j]. The initial conditions are given in the 
point (aci, y\,z\ ) reported in Eq. ( 5.15| ). 



not unreasonable that some of the nasty features of the present results might be changed if 
thick branes are included. In this case some of the unstable solutions might get stabilized 
since a new scale (the thickness of the brane) arises in the problem. For instance solutions 
with exponentially decreasing warp factors and decreasing dilaton coupling might turn into 
warped solutions with constant dilaton. On the contrary warped solutions with increasing 
dilaton coupling (corresponding to another class of fixed points emerging from the present 
analysis) might be driven towards naked singularities with Kasner-like features. 



VI. CONCLUDING REMARKS 

The effects of the dilaton and of higher order curvature terms are usually neglected 
in the context of six dimensional warped compactification. This is certainly a consistent 
assumption which can be, however, relaxed. In the present investigation six- dimensional 
warped compactification has been analyzed in the context of the gravity theory inspired by 
the low energy string effective action. Tree-level solutions have been derived. First order a' 
corrections have been also included precisely in the form induced by string amplitudes. 

The tree-level solutions exhibit a Kasner-like branch whose generic property is the occur- 
rence of curvature singularities. In this case the warp factors are powers of the bulk radius 
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p. The powers satisfy specific sum rules leading ultimately to divergences in the curvature 
invariants. The features of the tree-level solutions (as a function of the bulk radius) are the 
static analog of the time-dependent case which has been widely exploited in cosmological 
considerations. Kasner-like solutions can be found both for vanishing and non vanishing 
bulk cosmological constant. This means that, at tree-level, exponentially decreasing warp 
factors cannot be generically obtained. Moreover, since the Kasner-like behaviour leads to 
curvature singularities, there will be a regime where the solutions hit the string curvature 
scale where higher order string tension corrections cannot be neglected. 

The inclusion of the first order a' correction produces computable modifications in the 
/5-functions which can be viewed, for the present purpose, as a nonlinear dynamical system. 
Solutions with exponentially decreasing warp factors become then possible: they correspond 
to critical points of the dynamical system. Defining H = d p lncr and T = <9 p ln7 interesting 
critical points can be obtained in the case of constant (and negative) H and T . These 
critical points correspond to linear (decreasing) dilaton solutions. 

The obtained critical points are not always stable. It can happen that a given critical 
point correspond to stable/unstable node or to stable/unstable spiral points. In the case H = 
T the Kasner branch of the solution (for p ~ 0) can be analytically connected to a critical 
point. There are also solutions where two critical points can be analytically connected. In 
this case a given Tii < turns into 7i 2 ^ H\. Unfortunately in the numerical examples 
analyzed in the present investigation H2 > 0. 

In the case H 7^ T the situation is mathematically more complicated. By giving initial 
conditions of the system near a critical point singularities seem to be developed. Numerical 
evidence of this behavior has been presented by studying the singularity properties of the 
curvature invariants. 

It is now the moment of comparing what we have with what we ought to have. The 
ideal situation would be to find stable fixed points of the low energy (3 functions compatible 
with a constant dilaton field. In other words it would be nice to find solutions interpolating 
between a singularity and a critical point of constant dilaton or between an unstable critical 
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point and a constant dilaton solution. Such a behavior would give rise to exponentially 
decreasing warp factors with constant dilaton coupling which is exactly the original assump- 
tion tacitly adopted in five and six-dimensional warped compactifications. This situation is 
never realized in the set-up described in the present investigation but it cannot be excluded 
in principle in slightly different frameworks. It has been shown that exponentially decreasing 
warp factors are certainly compatible with the low energy string effective action, however, 
the obtained critical points do not have the wanted physical properties. These results should 
not be viewed as conclusive and might only be due to the limitations of the present analysis. 

We will now speculate about various possibilities which might cure the limitations of our 
approach allowing, hopefully, for stable critical points with frozen dilaton. As we discussed 
the first (obvious) limitation which should be relaxed concerns the dilaton dynamics. In 
the present paper only string tension corrections have been considered. In principle one 
can expect that corrections in the dilaton coupling should also be included for consistency 
since they become relevant roughly at the same scale at which string tension corrections are 
turned on. Following a complementary way of thinking one can also try to account for the 
effect of the dilaton potential which has been also (partially) neglected since we modeled 
it with a cosmological constant. An encouraging remark, in this direction, would be that 
already within our over-simplified set-up decreasing dilaton solutions were compatible with 
exponentially decreasing warp factors. These solutions might be analytically connected with 
a critical point of constant dilaton once a dilaton potential is present. A less encouraging 
evidence is however that constant dilaton solutions are obtained (within the present analysis) 
only as l/a! effect. These solutions are clearly non-perturbative and it might be that only 
by summing up all the string tension corrections a definite answer can be obtained. 

Another limitation of the present analysis is the absence of realistic (thick) brane sources. 
It is not unreasonable to expect that by adding a brane of finite thickness in the game some 
of the unstable critical points of the system will be stabilized at a constant value of the 
dilaton coupling fixed, presumably, by the tension of the brane. This possibility should be 
further scrutinized. Again an encouraging hint in this direction is that, already at the level 
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of our analysis, it is possible to connect a critical point to Kasner-like solutions. 

Furthermore in our analysis we completely neglected the possible presence of p-form 
fields whose effect should also be studied. To neglect the antisymmetric field strength is 
consistent if we look at the explicit form of the first string tension correction. However, 
an antisymmetric tensor field with large magnetic component might introduce qualitatively 
new solutions. In the present discussion the antisymmetric tensor field has been frozen. 
However, once an antisymmetric tensor is allowed on the same footing of the other light 
modes the structure of the string tension corrections changes radically: higher powers of the 
field strength will appear to first order in the string tension corrections. We leave all these 
themes for forthcoming investigations. 

Finally, some cosmological implications of our work should be mentioned. Once thick 
branes, antisymmetric tensor field and dilaton potetial will be taken into account in simplified 
models, time dependent solutions might be studied. In this case the problem will certainly 
be more complicated but not totally hopeless if stable critical points of the truncated (3 
functions will be found in the time-independent case. We also leave these topics for future 
investigations. 
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APPENDIX A: CURVATURE INVARIANTS 



In order to scrutinize the singularity properties of the six-dimensional metric discussed 
in the present investigation the quadratic curvature invariants should be properly discussed. 
The curvature invariants computed from Eq. ( |1 . 1|) can be written as 

h(p) = R^ aP R^ al3 = \{^~ 6ft 4 + AT 2 ? 1 + AT 12 + 16ft 2 ft' + 16ft' 2 ] (A.l) 
Hp) = C^C**"* = 1[2(W - T') + T{U - T)\ 2 (A.2) 

J 3 ( p ) =R = -[^JF 2 + m jZ + 1Q ^2 + 2jr / + m ^ ( A £j 

j 4 (p) = R^RT = l[F 4 + AHT Z + \m 2 T 2 + \m*T + 40ft 4 + AT 2 T' 

8 

+mTT' + m 2 T' + 4^ /2 + %F 2 H' + Sftj^ft' + 64ft 2 ft' + 167^'^' + 40ft /2 ] (A.4) 

The Weyl invariant vanishes in the case where a(p) and 7(p) are proportional, namely in 

the case where ft = T . The other invariants are do not vanish in the limit ft —>■ T . 

The curvature invariant are singular in the case of the tree-level solutions derived in Eqs. 

( |2.9D . For instance the Riemann and Weyl invariants can be written, for the solutions of 

Eqs. (F9|) , as: 

A 2 , 

Up) = — [4a 2 -12a 4 + /3 2 + 2/3 4 -4 (4a 3 + /3 3 ) cosh [v^ (p - p )] 
lb 

+ (4 a 2 + (3 2 ) cosh [2 V^A {p - Po)}} sinh [y/^L {p - p )f\ 

UP) = ^ (« - £) 2 A 2 (/3 - cosh [V=A (p - r )]) 2 sinh [</=A(p - p )f (A.5) 

This shows that the tree-level (Kasner-like) solutions are indeed singular for p — > p . Notice, 
again, that if a = (5 the solution still exists. In this case the Weyl invariant vanishes 
identically but the Riemann invariant (and the other invariants) are still singular. 

In order to detect singularities in a given solution the behavior of the curvature invariants 
has been scrutinized as a function of the bulk radius. This is a consistent procedure. In 
order to fully analyze the singularity properties of the solution we should also investigate the 
behavior of non-space-like geodesies (i.e. either time-like or null). This check is pleonastic 
if the curvature invariants diverge. However, the regularity of the curvature invariants is 
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a necessary but not sufficient condition in order to assess the regularity of a manifold. In 
principle, in order to claim that a given geometry is singularity-free we should chack that 
the curvature invariants are regular and that non-space-like geodesies are complete (i.e. 
they can be extended to any value of the affine parameter). These are the usual general 
relativity criteria for the analysis of physical singularities of a given space-time. Using 
Einstein equations together with topological considerations these criteria can be rephrased 
in terms of energy conditions involving the various components of the energy momentum 



tensor. Hawking-Penrose theorems are based on this type of considerations |31 . 

When singularities are analyzed in the context of the low-energy string effective action 
there is a further ambiguity related to the choice of the physical frame where the calculation 
is performed. In this paper the string frame parameterization of the action has been used. 
Other frames can be anyway defined. A very useful frame is the Einstein frame where 
the dilaton and the graviton are decoupled (at tree- level). Defining as a D = 4 + n 
dimensional metric in the string frame and as G^ v as the Einstein frame metric (in the same 
number of dimensions) we have that they are related through a conformal rescaling involving 
the string frame dilaton 0, namely: 

G>„ = e~^g^ u . (A.6) 

The (canonically normalized) Einstein frame dilaton (<£>) is simply related to <fi as <p = 
(f>\/2/ (n + 2). The transformation of Eq. ( |A.6| ) clearly induces a change in the curvature 
invariants. Consider, for instance, the scalar curvature in the string frame and call it R. In 
the Einstein frame the scalar curvature will be 1Z. The relation between R and 1Z is 



R = e 



-4> 



n-(n + 3)G a ?V a V^ - (n + 3) 4 (n + 2) G^V tt ^ , (A.7) 



where ip = ipy/2/(n + 2) and where the covariant derivatives are computed with respect 
to Gfj, u . The two frames certainly describe the same physics even if there has been debate 
on which of the two frames captures better the stringy features of a given solution of the 
low-energy /5-functions. As far as the singularity properties are concerned we can clearly 
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see that the two frames are not equivalent. Suppose, for instance, that the dilaton field is 
decreasing in the Einstein frame and suppose that IZ is regular in the same frame. Then 
Eq. (|A.7|) tells that R will be no longer regular in the string frame. Therefore singularities 
may appear and disappear from one frame to the other. 

This example only means that the two frames are indeed physically equivalent when we 
deal with physical solutions namely solutions where the dilaton is asymptotically constant, 
for instance. This observation helps in setting the goals of the present paper. We want 
warped compactifications, corresponding to stable critical points of the truncated /3-functions 
for which the dilaton field is frozen to a constant value. 



A similar situation |32| occurs in another type of inhomogeneous backgrounds (i.e. four- 
dimensional background geometries with two commuting Killing vectors which are hypersur- 
face orthogonal and orthogonal to each others). In |32] completely regular and geodesically 
complete solutions of the low-energy f3 functions have been discussed in the framework of 



these backgrounds. Also in the case []32| the crucial requirement is that the dilaton interpo- 
lates between two constant values. 



APPENDIX B: EQUATIONS OF MOTION WITH THE FIRST a' CORRECTION 

Inserting the metric given in Eq. (|1.1|) into Eq. (|3.4|) the reduced form of the action is 
obtained to first order in a', namely 

S = -J^ J d A x J dp d6 a 2 ^ e~+ {^- + 5H 2 + 2HT + 4d p H + d p T - (<9 P 0) 2 + A 

-e\m 2 T 2 + yft 4 + m 2 d p T+ YlHTdpH + XTH^T + 12H 2 d p H - (<9 P 0) 4 ]}. (B.l) 

Recalling now that 7i = d p lna and that T = <9 p ln7 we can perform, separately, the 
variation with respect to 0, a and 7. The variation with respect to gives 

-2 A - T 2 - 4 T H - 1 H 2 - 2 T' - 8 W + 2 T 0' + 8 H 0' - 2 /2 + 4 0" 

+e (6 T 2 H 2 + 24FH 3 + 15 ft 4 + 12 H 2 T' + 24 T H H' + 24 H 2 W 

-4 T 0' 3 - 16 U 0' 3 + 6 0' 4 - 24 0' 2 0") = 0, (B.2) 
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whereas the variation with respect to a and 7 gives 



fi{p) + 2(H+^- <t>')f'M - M ~ (W + f - 0U(p) 



+ 



(B.3) 



g'i(p) + 2(2^ - I - 0')^(p) - M ~ W ~ f - 0Xp) 



+ 



(2H- 



tA2 



(2W 



93(p)+9i(p) = 



(B.4) 



where the f(p) are given by 



t/4 



hip) = 2A-2FH+(J r + 2H) 2 + 2 (f' + 2H')-2f 2 + e (3 ft 4 + 12 ft 2 ft' + 2 < 
/ 2 (p) = 2 (.F+W)-6e [H 3 + ^H (2W + .F) +4HH' + 2 (HF + FH')] , 
/ 3 (p)=4[l-3eW(.F + W)], (B.5) 

and the g(p) are given by 



9i(p) 



1 

4 L 



2A + .F 2 -4.Fft + 10H 2 -2F + %H' -2<P' 2 , 
-e (6 ^ 2 H 2 - 2AJ r H 3 + 15 H A - Y2H 2 T' - 2AFHH' + 24 H 2 H' -2 /4 ' 
^ 2 (p) = _jr + 2 ^-6e (-FH 2 + 2H 3 + 2HH') , 

^ 3 (p) = l-6e^ 2 . (B.6) 



Recall that ' = <9 p . Eq. (|BT2|) exactly coincides with Eq. (|3T5|) . Inserting Eqs. (|BT5|)-(|B~6|) 
into Eqs. ( |B.3| ) and ( [B.2| ) we get, respectively, Eqs. ( |3.6| ) and (|3~T|). 



APPENDIX C: EIGENVALUES AROUND THE CRITICAL POINTS 



The two eigenvalues in the case of A < can be written as 



W 1 (k) 
W 2 {k) 



5A_ 



5A_ 



1 + 



212 k 



T(k) + y/S(k)Q(k) 



V(k) 

2Uk[ T(k)-^S(k)Q(k) 



V(k) 



(C.l) 
(C.2) 
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where 



1 r / 212 k ~\ 

T(k) = - 1 -1925 - 260180 k - 7570944 A; 2 + 5 W 1 + — ^— [385 + 384 k (-19 + 792 k)]j 

r / 212 k "i 

S(k) = 1 1015 + 99916 k + 2411712 k 2 - J 1 + — g— [1015 + 12 fc (2269 + 158400 k)} j 

r / 212 i 

Q{k) = 1 1 25 + 123380 A; + 5006592 A; 2 + 5 Wl + ^— [-25 + 48 A; (563 + 12672 k)} j 

V(k) = k [127925 + 576 k (12473 + 348480 k)] (C.3) 



Recall that k = eA_ 
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